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A transportation network is considered. The traffic demands associated with pairs of nodes and 
the (convex) traveling cost functions associated with the links are assumed given. The two problems 
of rinding the traffic patterns which either minimize the total cost <>r equilibrate the users" costs are 
formulated, and algorithms are constructed for the solution of these problems. 

Key words: Algorithm; least cost; traffic allocation; transportation. 

Introduction 

Many economic systems can be visualized as networks where nodes stand for commodities, 
and links and paths stand for simple and complex production processes. The type of system which 
can be thus described in the most natural way is probably a transportation network. In this case 
the nodes stand for "cities, " the links stand for roads directly connecting two cities, and the paths 
stand for roads connecting two cities directly or indirectly. 

A certain demand is associated with every pair of connected nodes of the network. This 
demand will be distributed among paths which join the pair of nodes. This gives rise to a traffic 
pattern, the determination of which is known as the traffic assignment problem. With every link 
of the network we associate a "traveling" cost which is assumed to be a function of the "traffic 
volume" on the link. We assume that the units traveling along this link uniformly share this cost. 

In some cases the traffic pattern can be regulated by some central authority, as for example, 
a network used for the transportation of military supplies or for a railroad network. It is obvious 
* that in this case, the problem which the central authority faces is to determine the traffic pattern 
which minimizes the total cost over the whole network. 

On the other hand a broad class of transportation networks can be described as user opti- 
mized. Here travel patterns are set up by individual users each choosing the cheapest way (in the 
light of other users' decisions) to arrive at his respective destination, rather than having his travel 
pattern dictated by a choice consistent with some aggregate system optimum. 

That the two above criteria lead generally to different traffic patterns was observed first by 
Pigou [1, p. 194] 3 in an example of a simple two node, two link network. Interest in this problem 
has been revived by Wardrop [2], who calculates the traffic patterns according to the above two 
criteria for the case of a network consisting of two nodes connected by n independent paths and 
for a special cost function. Wardrop discusses briefly the case of a general network and sketches 
the equilibrium equations, but he does not discuss their solution. 

Since 1952 several authors have reexamined the problem of flow patterns in a transportation 
network. For a complete bibliography we refer to a survey acticle by Beckmann [3]. We should 
observe here that two problems discussed by Wardrop, the problem of calculating the flow patterns 
according to the above two criteria, and the problem of planning an optimal investment allocation 
for improvement of the traffic network, still remain open. 
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Some progress towards the calculation of the flow patterns has been made. Almond [4] has 
constructed an algorithm for the solution of the user optimized network in the case of very simple 
networks. However, no extension of the algorithm for more complicated networks and no proof 
of convergence has been provided so far. 

A different method of attack is based on the observation that the user optimization problem 
can itself be reformulated as a total cost minimization problem for an appropriately chosen objective 
function [5, 3]. 

When viewed in this manner the problem is of the "multicommodity network flow" class, 
which has been considered [13, 14] in the literature. Tomlin [14] has shown that for the case of linear 
cost (congestion) functions, the problem reduces to a linear programming problem that can be solved 
fairly efficiently by the Dantzig-Wolfe decomposition principle. Others [11] have suggested the use 
of convex programming techniques to get around the nonlinearity of the objective function. In 
fact, it was the enormous number of constraints associated with the convex programming formula- 
tion of the problem for the simplest of networks that led us to develop the special algorithms pre- 
sented in the paper. 

Returning to the Tomlin algorithm, it should be emphasized that his algorithm takes advantage 
of the linearity of the objective function; on the contrary, the success of ours hinges on the non- 
linearity of the objective function, as will be demonstrated below. The algorithms should be viewed 
as a contribution to the theory of nonlinear multicommodity flow, as well as a contribution to the 
traffic flow literature. 

In the present paper we mainly try to solve the open problem of the calculation of the traffic 
pattern in a general network, for the two criteria proposed by Wardrop. Some progress has been 
made [10] on Wardrop's resource allocation problem; these results will be reported in a later 
publication. 

The paper is divided into two sections: the first concerns itself with problem formulation. 
Section 1.1 introduces the notation to be used, and the concept of a feasible flow pattern for a 
network. Section 1.2 describes two problems associated with transportation networks. The first, 
Pi, is to find a feasible flow pattern that minimizes the total cost of traveling in the network; the 
second, P 2 , to find the feasible flow pattern that would be arrived at if users considered only their 
own interests in choosing these paths. Section 1.3 spells out the conditions that are assumed con- 
cerning the congestion functions for the links of the network, and then goes on to give the necessary 
and sufficient conditions for the existence, uniqueness, and stability of a solution for problem Pi. 
Then the same conditions are derived for problem P 2 by showing that there is always a problem Pi 2 
associated with problem P 2 whose solution is that for P 2 , yet whose formulation is that of Pi. This 
theorem is a translation into the present paper's perspective and notation of the result of Jorgensen 
[5] referred to before, and implies that every traffic assignment problem of a user optimized network 
can be solved by solving the associated problem of a total optimized network. In addition, we 
extend Jorgensen's work by examining the stability of user optimized networks, as well as giving 
a more general condition for the user optimized and total cost optimized travel patterns to coin- 
cide. The section concludes with the conditions on the congestion functions that will cause the 
solutions of Pi and P 2 to coincide. 

Section 2 constructs algorithms for the solution of problems Pi and P 2 . Section 2.1 introduces 
the concept of an equilibration operator, and the conditions that must hold for such operators 
to obtain a solution to Pi, referred to as the process of "inducing an algorithm for the solution 
of Pi." The section concludes with the introduction of the notion of disjoint paths. In section 
2.2, we construct two equilibration operators, E ( t S j for networks with disjoint paths, and E„dsj, 
for any network. The operators are first applied to quadratic models, and we discuss under what 
conditions they induce algorithms for Pi. In brief, we show that Edsj induces an algorithm for 
simple (disjoint paths) and almost simple (see text) networks with quadratic cost functions, and give 
evidence that Edsj converges rapidly to a solution. Next, we show that the operator Endsj induces 
an algorithm for problem Pi for arbitrary networks with quadratic cost functions. 
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Section 2.3 extends the results of the previous section to cases where the cost function is 
required only to be twice continuously differentiable and convex, rather than quadratic. Section 
2.4 briefly compares the two operators and presents the respective conditions that appear favorable 
for their use. 

1. The Problem of the Traffic Distribution in a Transportation Network 

1.1. Generalities 

We start by introducing the concept of a transportation network. Let ^ be a network in the 
sense of Ford and Fulkerson [7, ch. 1, sec. 1], i.e., S? is a pair («#*, «=§?) where Jf is a collection of 
elements which will be called nodes and J? is a set of pairs of ordered elements of Jf which will 
be called links. 

By a path connecting the ordered pair w = (x 9 y) of nodes we mean a sequence of links (x\ , #2), 
(#2, ^3), • • ., {xn-u %n) where X\, x 2 , . . ., x n are distinct nodes, x\= x, and x n — y. Thus a path here 
is a chain in the terminology of [7, ch. 1, sec. 1]. In particular, every link is a path. The set of all 
paths of S? will be denoted by 8P. A pair w of nodes will be called connected if there exists at least 
one path connecting w. The set of all connected (ordered) pairs of nodes of ^ will be denoted by 2^. 
The set of all allowable travel paths which connect a w will be denoted by ^*. 

With every w = (x y y)eW we associate a nonnegative demand d w for travel with origin x and 
destination y. This demand will be distributed among all paths in 3P W . Suppose that petPw By 
fp we denote the part of d w which travels through /;. Thus we have the conservation equations 

(L= £./„. (1.1) 

We define 

A fixed value of ZF will be called a. flow pattern since it characterizes completely the flow. In the 
present paper we assume that the traffic flows are nonnegative real numbers and that the links 
of the network have infinite capacity. 

We will assume that a cost c (l is associated with every aeJ£ of S?. The value of c« is assumed 
to be a function of the total amount of traffic f„ through a. That is, 

Ca=C a (fa) (1.3) 

where 

/«= 2 W/> (1.4) 

with 

_ f 1 if a is contained in p, 
ap \ otherwise. 

We define 

W= {f a : ae^} y % = {c a (fa) : aeJ?, / a€ [0, 00)}. 

The triple ST = {& y 2), &} will be called a transportation network. 

Throughout the paper we consider problems of the following type: A transportation network 
3~ is given and the flow pattern J^ is the basic unknown. So f ar <jF has to conform to the conserva- 



4 We assume that ^ u . is a reasonably small set which can be enumerated in advance with little difficulty. This assumption, certainly a plausible one to make foi 
traffic networks, avoids the problem of computing all paths in a network, an enormously time consuming task for large networks. 
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tion equation (1.1). An 3F which satisfies (1.1) will be called a feasible flow pattern. The set of all 
feasible flow patterns (for fixed S?, 2) will be denoted by <2?[S^, 2\ It is obvious that there exists 
a unique feasible flow pattern only in the case for which for any weW which is connected by more 
than one path, d w = holds. Leaving aside this trivial case, we observe that there is an infinity of 
feasible flow patterns. 

We are now ready to formulate the two basic problems with which we will deal in the paper. 



1.2. Formulation of the Problems Pi and P 2 

PROBLEM Pi [5^]: Given a transportation network ZT = (^, 2f y ft), find a feasible solution 
1F\(?F) which minimizes the total cost 

C(F) = 2c„(/L) (1.6) 

spent in the network. 

A solution 1F\(?F) of problem P\\ST\ will be called a "system optimizing" flow pattern. 

As noted in the introduction, this is a reasonable problem but in many cases the network is 
in fact "user optimized." Each user of a link a will be charged with a portion of the total cost c n 
on this link. It is natural to assume that there is full interaction between all units traveling on 
link a; that is, the cost is distributed uniformly among them. Thus, the share of the cost of each 
unit traveling on a will be given by 

<?«=««(&)- ^ (L7) 

fa 

In consequence, the personal cost c p of a unit traveling on pe& will be given by 

Cp = ^8apC a (1.8) 



where the incidence symbols 8 ap have been introduced by (1.5). 

In order to make clear the notion of a flow pattern which is "user optimized," we introduce 
the following definition. 

DEFINITION (1.1): For given 2T= (fg, 2, %>), by an equilibrium flow pattern IF' we mean a 
feasible flow pattern with the following property. Let weW such that d w > 0. Choose any pe8? w 
for which fp > 0, and any number A/, < Af<f' p . Consider another path qe^ w . Then the individual 
cost c p (W)Afof A/ in the original flow pattern 9>* is not greater than the individual cost c q (W r ') A/ 
in the flow pattern <jF" defined by 

/y=/Z+A/, (i.9) 

f;=f' r , , re^, r^p, q. 5 

In other words, an equilibrium flow pattern is an equilibrium point in the sense of Nash 
(e.g. [8, sec. 7.8]) of the noncooperative game among the various users of the network. Having 
given the definition of an equilibrium flow pattern, we now formulate problem P 2 . 

PROBLEM P? [&~\: Criven a transportation network «^~, find an equilibrium flow pattern 

1^2 = ^2 OT). 



5 Under the assumption 4? to be imposed shortly, this condition implies that no distribution of A/ among several paths of ^ u .— {p} reduces the cost for A/. 
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1 .3. Study of the Solutions to Problems Pi and P> 

It is not to be expected that the problems /\ P 2 , formulated above are well posed unless some 
conditions are set on the form of the cost functions c a (fa)- Whenever we consider the problem l\ 
we will assume that the above functions satisfy the following assumptions for all aeJ£ . 

1. c a (fa) is continuous on [0, oc). 

2. c a (0)=0. 

3. c a (fa) is strictly increasing on [0, °c). 
4i. c (l (fa) is strictly convex on [0, oo). 

The interpretation of conditions 1—3 is obvious. Conditions 4i, for differentiable c a (fa), means 
that the rate of increase of the cost, i.e., the marginal cost, is a strictly increasing function of the 
traffic flow f a (congestion effect). 

Whenever we consider the problem P 2 we will assume that conditions 1—3 above are satisfied 
but, in the place of 4i, we will impose the condition 

4-2. c a (fa) is strictly increasing on [0, o©) 

with the understanding that 

c„(0)= lim^- 
/;-><>' fa 

Assumption 4 2 provides a slightly different interpretation of the congestion effect with the emphasis 
placed on the individual rather than on the marginal cost. In fact conditions 1, 2, 3„4| imply condition 
4 2 . (However, conditions 1, 2, 3, 4 2 do not imply, in general, condition 4|.) 

The simplest model which satisfies the above requirements corresponds to a cost function of 
the form 

Caifa) =gaM + hafa, go > 0, h a ^ (1.10) 

and will be called the quadratic model. In this model the congestion effect depends linearly on the 
traffic flow. 

Having specified the admissible form of ^, let us consider the problem l\ [< < T],J r = (^, 3 y <€). 
Recall that this problem calls for the vector (<-F, ~W) which solves the minimization problem: 

minC(W)=y / Ca(f a ) (1.11) 

subject to 

pW 

/ p S0;p€^, (1.12) 

pe&w 

Observe that a fixed ^ induces a unique 3F through (1.12)i. But it is possible that more than one 
feasible J^ induces the same 3F . The set of all feasible 3P which induce a given fixed !F will be 
denoted by R[3^]. On account of (1.12), R[^] is a convex set. ^ will be called feasible if R\_.W\ is 
non-empty. The set of all feasible ^ will be denoted by JT[^, £?]. In appendix I we give an example 
of a transportation network such thai for someFeJ[^, ^], R[^F] contains an infinite number of 
elements. 

Note that if (3F\ 3P'), (^", ^") satisfy the constraints (1.12) so does any convex combination 

(&*¥) = \'(^\W) + X"(:?\ ¥"), k\ \">0, \' + A"=l. 
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On the other hand and on account of 4i , 

C{¥) ^k'C{¥') + \"C(¥") 

and equality may hold only if W* = W\ Consequently, P\ is a convex minimization problem in 
3£®3£ and, in particular, it is strictly convex in 2T. Using the theory of convex programming and 
observing that the total cost function depends only on ~W , we arrive at the following theorem. 

Theorem (1.1): Given ST = &, 3), <g), there exists a unique tFieW^g, @j\ such that C(3P\) is the 
minimum of C{?F) over 2£\3 , Qi\ Every element JSe/?[<jFi] is a solution of problem P\. 

Thus, problem P\ always possesses solutions and in particular it possesses a unique solution 
if and only if R[^\] consists of a unique element. 

In the special case where the c (i {f (l ) are differentiable functions we can prove the following 
theorem. 

THEOREM (1.2): The flow pattern ZFe3£ is a solution of problem P x if and only if it has the 
following property. For any weW connected by precisely the paths p u . . ., p m , these paths can 
be so numbered that 

4,(F)= . . . =c; i .(F)-M u .g c ;, s+i (F)g . . . sc;jW), 

f Pr >0, r=\, . . ., 5, (1.13) 

//< =0, r=s+l, . . ., m, 



where we use the notation 



c;(^)-2s„ /)C ;,(/„), 

ae.2? 

dCa{fa) 



C' a (fa)=- 



dfa 



PROOF OF SUFFICIENCY: Assume that ^e£ satisfies (1.13). Let ^+A<FeiT be a feasible 
reallocation. The change of the total cost is given by 

AC=2fc»(A + A/,)-c,(/«)]. 

Applying the mean value theorem and using the fact that the functions c' a {f (l ) are (strictly) increasing 
we obtain 

ACs^cX/JAAr A^ (1.14) 



Recalling (1.4), 



Then, 



A/-„ = 2 8 «"4/p. (1-15) 



AC g £ 2 8apC' a (faWp= £ ^fpCp(^). (1.16) 

Note that if/ p =0, then A/ p ^ 0. Thus, using (1.13), 

Af P c' p (,W)^Af p M w (1.17) 
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where w is the pair of nodes which is connected by p. Recalling that y A /,, = (). W( ' obtain from 
(1.16) and (1.17), /f/ " 

AC i= 

which proves that J^ is a solution of problem f\. 

Q.E.D. 

Proof of NECESSITY: Suppose that J^e^T is a solution of problem Pi, but there exist paths 

p, <7€^V such that/p > and 

c' p (<F)-c' q (W) = e>0. (1.18) 

Assume now that a portion A/of/p is reallocated to the path q. The change of the total cost is given 

by 

AC = £ 8UCa(fa " A/) - C a (fa) ] + £ 8g„ [c a (/« + A/) - C a ifa) I 

a€ jgf a€.S? 

where 

s | Sap if a is not contained in q, 



'<!/> 



if a is contained in q, 



and 8#, ; is defined in an analogous fashion. 

Applying the mean value theorem and recalling that c' a (fa) are (strictly) increasing functions 
we end up with , 

AC < j - 2 8fcci(£ - A./) + ^ 8g,c;c/L + A./) a./: 



Now V 8gpC a (/ a — A/), y ^^/h + A/) are continuous functions of A/. Hence, we may choose 
a positive A/ (for feasibility it must be such that /,, — A/ ^ 0, whence A/^/p) such that 



Hence, 



2 s,Vv',</«- A /)> 2 W/«)-7. 
2 %,r,;(/;,+A/)< 2 8 /«/ c «</«>+ € 3 

AC < - £ 8fa' a (faW+ \ A/+ 2 o^(/«)A/+f A/. 
It is easily seen that the above inequality may be written in the form 

AC < f - 2 SapC'aifaH £ «**(£) +Tf } A/' 



or. using (1.18), 



ac<-|a/<o 



which is a contradiction to the assumption that J^ is a solution of problem Pi. Q.E.D. 

Actually it can be shown that (1.13) are simply the Kuhn-Tucker conditions (see [9, ch. 6]) 
for the minimization problem (1.11), (1.12). However, these conditions have been derived inde- 
pendently here in order to keep the paper self-contained. 

In the case of the quadratic model. P\ reduces to a quadratic programming problem and 
(1.13) become linear. 

From the convexity of problem P\ we can obtain additional information, namely that the 
solution is stable. To make this precise we introduce the following definition. 
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DEFINITION (1.2): Let e > 0. We say that ;¥e^E satisfies (1.13) modulo e if and only if for any 
p, q connecting any weW and such that f p ^ de, either 

\c\fW)-c' q (¥)\^- d e (1.20) 



holds, or else both 

d 



c' v (^)<c' q (^)+-e (1.21) 



and 

f q < de 

where c, <i are arbitrary but fixed magnitudes having dimensions of cost and traffic flow, respec- 
tively, and included in order to make e dimensionless. 

Theorem (1.3): Let & ' — (&, Q)^) be a transportation network with twice continuously differ- 
entiable cost functions, and ,¥ ' a solution of problem Pi[^"\. Then there exist numbers K and L 
which depend solely on 3~ , such that 

C(W)-C(W)< cKe, (1.22) 

\\W-¥'ft - 2 \f«-fa\ 2 < d 2 Le, (1.23) 

for any 3Fe3£\2T\ which satisfies (1.13) modulo e, e > 0. 

PROOF: Assume that J^eJT satisfies (1.13) modulo e. We set A^ = ^' — 3F . As in the proof 
of Theorem (1.2) (see eq (1.16)), we have 

C(W f )-C(W)^^Af p c;(F). (1.24) 

p€& 

We decompose 3T,3r = 3TiU3r 2 . so that if weW u then f p < de for all pe^ w , while if weW 2 , 
then f }) ^ de for at least one pe!P w . In particular one of W\, W* may be empty. 

For weW\ 

Af p >-de,pe0> u , 

Hence 

^^^(Fl^-c-eJS^Fl^-c/de (1.25) 

where K\ can be easily estimated in terms of elements of ST . 
Fix now weW 2 , and let ^ w = {p u • . •, P™ u }- Suppose that 

f Pr ^ de, r= 1, . . ., s, and f Pr < de, r=s+ 1, . . ., m w . 
We write 

l*p r = Cp r (^) — c Pl (^),r=l m, r 

and we observe that since J^~ satisfies (1.13) modulo e, 

i , c , 

\Pp r \ <-€ , r= 1, . . ., 5, 
a 

(JLp> e ,r=s~\-l, . . ., m w . 

d 

It is also easy to obtain an estimate of the form 

fXp^-E w , r=5-hl, . . ., rciw 
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where the number E w depends at most on J and w. For example 

■ £^max r J] 8 aPr C 'a{ S *° 



c 

d 



(1.26) 



where // „ is the set of all weft •> which are connected by at least one path containing a. Notice also 
the useful estimates 



He 



Af„ r > — <ie ,r = s+\.. . ., m w , 
|A/ Pf .| S</„ , r=l, . . ., m ir . 

£ A/pC p (^)=c Pi (^) 2 A/.+ S^,A/„ =2 ^ r A/ Pr + X M» P 4&, 

/*./„■ />€./„■ / = i r=l r=s+1 

^ | — sc-j— (m w — s)c max 



<*«, 



('^AV jl^-^niax^, 



From this last inequality, (1.25), and (1.24) we deduce (1.22) for 



d w 



K = ^ /;/ " max ( "7/ » ^ ,r ) "*" ^ , - 



(1.27) 



We now proceed to the proof of (1.23). Note thai 

C(W) - c(W' ) = X ^7 A/;, + ~ 2 fl "7T? ^-^" 

with — = calculated at the point &' and — = — -calculated at a fixed intermediate point &' + 0AJ*\ 
d/a dfadfb 

^ ^ 1. Since #"' is a solution of problem P, and ^' + A^ r =^e&, then 



ac 



a/;,=o. 



On the other hand 



Hence, 



d^ ( c a(/«), if a =6, 



d/«d/& (o , if a ^ 6. 



C(^)-C(&")=-^cZ(fa)btf' 

tie y 

Since c„ (f a ) are strictly convex, the constant 

— k = min min c'a(f a +0Afa) 

U~ aeJ/ 0e|(). 1| 

99 



(1.28) 



is positive and (1.23) follows with the help of (1.22) for 

L^~ (1.29) 

Q.E.D. 

The interpretation of the above theorem is that the solution of problem Pi is stable. We em- 
phasize here that, as follows from the proof, the constants K, L can be estimated explicitly in 
terms of known characteristics of the transportation network. Then, apart from its theoretical 
importance, the Theorem (1.3) is useful in practice since it provides a means of estimating the 
distance of a given feasible flow pattern from the solution of problem Pi. An explicit application 
of the above observation will be presented in part 2 of the paper. With these comments we complete 
the study of problem Pi. 

We now proceed to a similar study for problem P 2 . We start by proving a theorem analogous 
to Theorem (1.2). 

THEOREM (1.4): The flow pattern J^eiF is a solution of problem P 2 if and only if it has the 
following property. For any weW connected by precisely the paths p\ 9 . . ., p m , these paths can 
be so numbered that 

c Pi m= . . . =c P8 (W)=A w ^c Pg+i (W)^ . . . ^c p jw) r 

/ Pr >0, r=l, . . ., 5, 

/p r =0, r=5+l . . ., m. (1.30) 

PROOF OF SUFFICIENCY: Assume that J^€ JT satisfies (1.30). Let p, q be two paths connecting 
the same we// and such that f p > 0. By (1.30) 

c p (W) ^c q (W). (1.31) 

Suppose that a portion A/, of f p , < A/ =f p , selects the path q. By 5^ ' we denote the resulting 
flow pattern. Recalling (1.19) we have 

C q (W')-C q (W) = 2 Eqitaifa + M-eaifa)] > 

where use has been made of the fact that c (l (f a ) is a strictly increasing function. In particular, 
recalling (1.31), 

c q (W) >c p (W) 

which shows that J^ is an equilibrium point in the sense of Definition (1.1). 

The proof of necessity is essentially a repetition of the proof of the necessity in Theorem 
(1.2) and will be omitted. , Q.E.D. 

The above conditions are nothing more than the standard average cost equality conditions 
an economist would expect to find in a system that is optimized by individuals acting independently 
of one another with no regard for total system optimization. Conditions (1.30) are known (but not 
in full generality) at least since the time of Pigou's treatise referred to in the introduction. Many 
authors consider the conditions themselves as a definition of problem P 2 . 

Comparing (1.13) with (1.30), we observe that there exists a remarkable similarity between 
them. The role of the average cost c p in (1.30) is played in (1.13) by the marginal cost c p . Starting 
from the above observation, we will now show that there exists a close relationship between the 
set of problems Pi and the set of problems P>. 
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Definition (1.3): Given P 2 [^] = P 2 f^, '/, 2 #], we construct a set of cost functions 2X % in the 
following way. For any aeJ? we set 

2lCa(fa)= \ a £alf)df. (1.32) 

Jo 

Note that if 2 c a satisfies conditions 1, 2, 3, 42, then 2X c a satisfies conditions 1, 2, 3, 4i. Furthermore, 
2\c (l is continuously differ entiable. The problem P\\^, 3)^ 2 \€\ will be called problem P x associated 
with problem P 2 [&> 2 ', &] and will be denoted by P\ 2 \fS y 2$, >&]. 

Similarly, given P x \fS, 2), \€\ where x ^ consists of continuously differentiate functions , we 
construct the set of cost functions X2 € through the use of 

12Ca(fa)^lC' a (fa)fa. (1.33) 

Note that if iC a (fa) satisfies conditions 1, 2, 3, 4i, then \ 2 c (l {f (l ) satisfies conditions 1, 2, 3, 4 2. The 
problem P 2 \f§ , £3 , 12^] will be called problem P 2 associated with problem Pi \fS , 2), &] and will be 
denoted by P 2X [&, 3, 1%]. 

The above given definition of the associated problem is justified by the following theorem. 

Theorem (1.5): Let .W be a solution of problem P 2 [^, 2 y <€\. Then ZF is also a solution of 
problem P V2 [^ y 0, #]. Similarly if .W is a solution of problem Pi[&> 2, &] with <g consisting of 
continuously differentiable functions , then ^ is also a solution of problem P 2X \ C 3 , r Z, ^|. 

PROOF: The proof follows from the construction of the associated cost functions. In fact we 
observe that (1.30) written for P 2 [^, ( /\ i&\ and (1.13) written for P x \^ y 2, 2 i^| coincide. Similarly 
(1.13) written for P x \<0, £d, x2 <g\ coincides with (1.30) written forP 2 [^,S, i 2 #|. Q.E.D. 

The notion of the associated problem is very simple but it will be of essential importance 
throughout the paper. For example note that Theorem (1.1) and Theorem (1.3) immediately imply 
the following corresponding theorem for problem P 2 . 

THEOREM (1.6): Given &' = (&> 2, 2 #) there exists a unique tF = ^F 2 (&, 2, >{& ) such that every 
^eR[3^] is an equilibrium solution of problem P 2 \fS^3l y2 €\. Furthermore, this solution is stable 
(in a sense quite analogous to the notion of a stable solution of problem P x induced by Theorem (1.3)). 

PROOF: Consider the problem P X2 [<& y 2, <$ ], i.e. the problem Pi [^ 3 % i 2 #] , and let F = Fi 
&y &, 12^) be its solution. Obviously, ^ is also the unique and stable solution to problem P 2 [&, 

OE.D. 

We will close this section with certain simple observations about the associated problems. 
In the case of a quadratic model, the associated problem is also quadratic. More precisely, if 
iCa(fa)=igcfa + hcf a then i 2 c a (/o) = g(f 2 a+ h t J„. Similarly if 2 c a (fa) =g(J? + h ( f„ then 2 \c (l (f (l ) 
— kgcfa + hqfa- In general, it is obvious that 212^ = 2^ and 121^—1^. 

A natural problem is the following: Suppose that a network ^ is given. Characterize the type 
of ^ for which the solutions of problem P x [^, 3) y %?] coincide with the solutions of problem P 2 [^, 
Sly tf] for every 3. Such cases are extremely desirable because in them the pattern created by 
the individuals acting in their own self interests coincides with the pattern most economical for 
the total society. We have already seen that the solutions of the associated problem P2[^> 3)<> &] 
coincide with the solutions of the associated problem P\ 2 \& y 2 y ( &~\. In consequence, the solutions 
of problem P x [&, 2,%?] will coincide with the solutions of problem P 2 [^, 2, &] if and only if they 
coincide with the solutions of problem P X2 [^, 2, &]. 

Recalling (1.32), we conclude that the solutions of problem P\[&, 3> &] coincide with the 
solutions of problem Pi 2 [^, 2, &] if 



nh _ _ 
c a (fa) = V ca(f)df. aeJ? 
Jo 



(1.34) 
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where rj is an arbitrary positive constant, the same for all aeJif. Integral equation (1.34) has the 
following solution: 

Ca(fa)=C%n (1.35) 

where c% is an arbitrary constant. In order for c a (f a ) , as given by (1.35), to satisfy conditions 1, 2, 
3, 4i or 4 2 . we restrict c%, tj so that c ( J ; > 0. r/eJ?f , n > 1. Actually (1.35) gives the most general type 
of cost functions which guarantee coincidence of the solutions of Pi\fS, 3, %?] with those of P 2 [^> 
&> %?] for an arbitrary ^ . It should be noted, though, that for special networks the class of such 
functions can be broadened. 

2. Development of Algorithms for the Solution of the Problems Pi and Pi 

2.1. Generalities 

In this section we develop algorithms for the solution of problem Pi[&*~\. Obviously if such an 
algorithm is available, the solution of problem Pi\ST\ can be obtained as the solution of a Pi prob- 
lem, namely the associated problem p2i[*^"]. 

Roughly, the method of solution can be described as follows: Starting from an initial feasible 
flow pattern we construct a sequence of feasible flow patterns which converges to the optimal 
solution. 

To be precise, we introduce the notion of an "equilibration operator." 

A map 

E w :&^>&, weW 

will be called an equilibration operator associated with weW if it sends ^e££ into <jF'eiF such that 

n=f„ 



unless pe&w 
A map 



£:^->iT 



will be called an equilibration operator associated with a transportation network <fF, if E can be 

factored, 

E = E Wn . . . 0E Wi (2.1) 

where {w\, . . ., w„} = W and E Wi is an equilibration operator associated with w\. 

We now give our definition of an algorithm. 

DEFINITION (2.1): Let 2T be a transportation network and E an equilibration operator associ- 
ated with 2T . We will say that E induces an algorithm for the solution of problem Pi \J?"\ if for any 

yW-*& 1 \JDi n^oo (2.2) 

where 

.fWspjTW n=l,2, . . . (2.3) 

and !W\ is a solution of problem Pi. 

The following theorem gives sufficient conditions for an equilibration operator to induce an 
algorithm for the solution of problem P x . 

THEOREM h (2.1): Let 9 be a transportation network and E an equilibration operator associated 
with 3~ and having the following properties: 

(1) EZF = !F for some ^eS implies that ^ satisfies (1.13) for all weW, so that W=W V 

(2) E is a continuous mapping from S to 2£ . 

(3) C(gg) ^ C(W)for all ^e$?. 

(4) C(E^)=C(^) for some ^e^ implies that E& ' = & \ 



s A proof of convergence along similar lines has been communicated independently to us by W. A. Horn of the National Bureau of Standards. 
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Then E induces an algorithm for the solution of problem P\[&~]. 

PROOF: Let ,^ (()) e^ and J^ ( "> = E",9> {{)) , n= 1, 2 We have to prove that 

W {n) ^Wi(,T), rc-»oo. (2.4) 

We first prove that every convergent subsequence {«^ n *>} of {^ {n) } converges to a solution J^i 
of P, [,r]. In fact, lei 

&r(n k )->&r t &-»oo. (2.5) 

Since 3£ is closed, J^ei?'. The sequence {C(J^ (;/) )} is decreasing and bounded from below by 0, 
hence is convergent. By Cauchy's theorem, given e > 0, 

|C(eF ( »a-)) - c(F { »k)) | = |c(3f ( "* +1 >) - c(W<*k))\ < | 
if Aa* L (€). 

Since C(J^) is continuous, lim C (&(**)) =€(&). Hence 

|C(<f< w *>)-c(F)|<| 

if A ^^(€). Furthermore, the continuity ofC(J^) implies the existence of 8 t some such that 

\c0w)-c0w')\ <|, if |F-^'| <s 6 . 

On the other hand, since E is a continuous mapping, given 8 > there exists 17(8) such that 

IfiF-H^'l <S i£\W-W'\ <y(8). 
Finally, from (2.5) it follows that given 77 > 0, there exists kn(rj) such that 

|^_^u A j| <T ^ lfk^k 3 (7)). 

Suppose now that k ^ max {Ai(e). //2(e) . A : ? ( 17 ( 8 e ) ) }. Combining the above results we obtain 

\c(eW)-C(W)\ ^ |C(£^)-C(z^< n *>)| + |C(I^ <e. 

But the left-hand side of the above inequality is independent of k and hence 

C{E^)-C(¥)=(), 

whence 3^ is a solution of problems Pi by (4) and (1). 

We now proceed to the proof of (2.4). Suppose that it is false. Then there exists a positive 
number 8 and a subsequence {^ (u k ) } such that 



1 1| = 



2IL/h' , -/>.,l>8. (2.6) 



The sequence {J^"/,*} is bounded. By the theorem of Bolzano-Weierstrass there exists a con- 
verging subsequence {^ ( "A/ } }. As proved above 

$frin kl )^>$r u /_> x 
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where 3F \ is a solution of problem P\. In particular, 

and this is a contradiction to (2.6). 

Q.E.D. 

The above theorem provides a criterion for establishing that a given equilibration operator 
induces an algorithm for the solution of problem P\. A limitation of the usefulness of the theorem 
may arise from the fact that it is not always easy to check whether assumptions 1-4 are satisfied. 
The following proposition simplifies this problem. 

THEOREM (2.2): Let {E w : weW} be a collection of equilibration operators associated with the 
pairs of connected nodes of a transportation network & '. Suppose that for every weW^ E w satisfies 
the following conditions. 

(1') E w & ' — £F for some J^eST implies that 3F satisfies {1.13) for this fixed w. 

(2') Ew is a continuous mapping from ££ to 2£ . 

(3') C(E^_)^C(W) for all ^e^. 

(4') C(E^W) = C(W)for some &e% implies that E w & = &. 
Then any equilibration operator associated with ZT and constructed by composition of the above 
collection {Et w :weW} satisfies conditions 1-4 of Theorem (2.1). 

PROOF: Assumption 1 follows easily from 1' and the structure of an equilibration operator 
associated with a pair of nodes. Assumption 2 is an obvious consequence of 2'. Similarly 3 follows 
immediately from 3'. Finally 4 follows by a combination of 3' and 4'. 

Q.E.D. 

The above theorem reduces the problem of checking conditions 1-4 of Theorem (2.1) to the 
much simpler problem of checking conditions 1'— 4' of Theorem (2.2). 

Sometimes an equilibration operator E associated with a transportation network satisfies 
conditions 1 and 2 of Theorem (2.1) but it does not satisfy (or at least we cannot prove that it satis- 
fies) conditions 3 and 4. Then of course we do not know whether E induces an algorithm for the 
solution of Pi[&~], Nevertheless the sequence {E"^ {0) } may lead to the solution of problem 
P\\JT\ as shown by the following theorem, the proof of which is similar to the proof of Theorem 
(2.1). 

THEOREM (2.3): Suppose that an equilibration operator E satisfies conditions 7, 2 of Theorem 
(2.1). Suppose further that for some choice of J^ (0) the sequence {E "#" (0) } converges as fi'—»o°. Then 
{E n ^ {i)) } converges to the solution ^y, of the problem P x . 

REMARK (2.1): We have seen that an equilibration operator E which induces an algorithm for 
the solution of problem P\ enables us to calculate through (2.2) the unique $F\ associated with a 
problem P\[&~]. Then we know that R[3^\] is the set of solutions of problem Pi. The calculation 
of an element of /?[J^i], given J^i, amounts to finding a solution to the system (1.1), (1.4), which 
might be accomplished by phase 1 of the Simplex method. This requires a rather tedious calculation. 
However, as shown in the proof of Theorem (2.1), some elements of /?[J^i] can be obtained directly 
from the algorithm as limits of the convergent subsequences of {£' / 'J^ (0) }. In particular, if R[^i] 
consists of a unique element then 

E n 9 rw *-*9 r u n^>^- 

REMARK (2.2): The stability results of Theorem (1.3) can be employed here in order to estimate 
| \^\ — E n ^ (0) \ | , n=l, . . . , and thus they provide a means for deciding whether the approxi- 
mation is satisfactory, in which case the algorithm is terminated. In fact the estimation of the 
smallest e modulo which E n ^ m satisfies (1.13) can be obtained for example by a method com- 
municated to us by Alan Goldman and which is given in appendix 3. Then a use of estimates 
(1.22), (1.23) reveals the accuracy of the aprpoximation. 

Recall that the proof of Theorem (1.3) involved an estimate of the form 
2f]2 

||£,^(o)_j^ f^~ [C(£»^ (0) )-aF,)], 
104 



which enables us to estimate the convergence of {E n ,W (Q) } in terms of the convergence of 
{C(E",^ {{}) ) }. Thus when we apply the algorithm it is sufficient to inspect the sequence of the 
successive values of the total cost. From the rate of convergence of this sequence we can judge 
the rate of convergence of the sequence {E n ,j^ ({)) }. 

It is convenient to introduce here the concept of disjoint paths. Let weW. The set & w will be 
called disjoint if there is no ae'/ which is contained in more than one pe^, r . 

A network .0 will be called simple if tP w is disjoint for every we^ w . 

In the following paragraphs we will construct two equilibration operators Edsj, E n d.y and we 
will discuss under what conditions they induce algorithms for the solution of problem /V These 
operators will be introduced first for the quadratic model and then the definition will be extended 
to the case of general convex cost functions. In particular Edsj can be applied more naturally to 
simple networks while Ends} has been designed for application to nonsimple networks, for which 
Edsj need not induce an algorithm for the solution of problem P\. 

2.2. The Quadratic Model 

Let ,f — {//^ /, 9o) be a transportation network with quadratic cost functions. 

c„(f a )=gji+hj„ 9 aeJf. 

According to (2.1), in order to define the equilibration operators Edsj, Endsj associated with the 
given ./ it is sufficient to define their factors E$&, E'," is J. respectively, for every weW\ 

a. The Equilibration Operator Ed$j for the Quadratic Model 

We start by motivating the definition of EftK Let /ce// such that ,/V [p u . . .*/>,„} is disjoint, 
and consider any ./ 7 e/. By %$ we denote the subset of 2 such that 

Let us seek the element ,f e / ,; which 

minimizes C(,F' ) over ) w ■ (2.7) 

In order to solve this minimization problem it is convenient to introduce the following notation. 

£ r / V = X 8 <*Pr#«> (2.8) 

at / 

M/ V (^> s 2 8ap r [2gJ a + ha] , (2.9) 

hp r (^)^fM Pr -2g p / Prt ( 2.10) 

Without loss of generality assume that 

//„,(./' )^ . . . ^h p {&). 

m 

A comparison with (1.13) leads to the conclusion that the solution J^"' of (2.7) satisfies the 
following equilibrium condition: 

f-i IJ} + 2g p , (jp\ —f p j) = . . . = /x ; , s + 2g Ps (f Ps —fp 8 ) = l*>p s+x + 2gp s + 1 ifps+i ~fp s +\ ) 

= • • • = l*<P m + 2gp m ifpm ~fl> m ) » 

£ r >0, r=l, . . .,s ;j £ r =0,r=s+l, . . .,/». (2.11) 
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On account of (2.10), (2.11)i reads 

M w = 2g Pl f Pl + h Pl = . . • = 2gp s fp s + h P8 ^2g PM f^ 1 + h P8+1 ^. . . ^2g Pm f Pm + h Pm 

and the solution of (2.11) gives 

M w -h„ r 
fp ~~2~i ,r=1, • • -' s ' 



fp r =Q, r=5 + l, . . ., ra, 



yhere 



(2.12) 



2d w +^h p Jg Pr 

M w = ^ . (2.13) 

r=l 

Thus «!F' may be calculated through (2.12), (2.13), provided that the critical index 5 is known. 
We now give a procedure for the evaluation of 5. From (2.11), (2.12) we obtain the condition 

h Pm ^ i= h Ps+1 ^ M w > h P8 ^ h P8 _, ^ i= h Pl . (2.14) 

Let 

f 

2d w +Yhi>Jgp k 
M' w = ^ ,r=l, . . .,m. (2.15) 

A = l 

The index 5 for which Mf r satisfies (2.14) is the critical one. The existence of a unique s having this 
property follows from the existence ol a unique solution to the minimization problem in question. 
Nevertheless, we demonstrate separately the existence and uniqueness of such an s, using the 
following identities which are of interest in themselves: 

Mi=YJ^-+h Pl >h Pt , (2.16) 

{Ml- /W'„r') § Ug Pk = llg Pr (h Pr -M r w ) , (2.17) 

A-=l 

or 

(M^-M^^llg.^llg^h^-M^). (2.18) 

Let S be the set of indices such that re S if and only if M r w > h Pr From (2.16) it follows that 
US. Suppose that reS. Then M r w > h Pr . From (2.17) M^r 1 >M r w . Thus h Pr _ l ^ h Pr < M' w < M r u r\ 
i.e., (r— l)e S and hence 1, 2, . . ., reS. Let s be the maximum index in S. Then h Pg < M s w . Further- 
more, either s = m or s + l^S which implies h Ps+l ^ M** 1 . Using (2.17), (2.18), we conclude that 
M^ 1 ^ h Ps+l implies M 8 W = M s ^\ which in turn implies M s w = h p v Thus the existence of a unique 
s has been established and another method of construction of 5 (as the maximum index in S) has 
been given. 

Summarizing, to calculate the solution ,F ' of the minimization problem (2.7) we apply the 
following procedure: 

(1) We calculate the quantities g p , h p (J^), r= 1, . . .,/?/. 

(2) We arrange h p [fW) in nondescending order and we relabel them according to this order. 
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(3) We calculate M r w , r=l, . . ., m from (2.15). 

(4) If M , £y^hp m we set s = m. If M$ ^ A Pm , we find the unique index s such that h l)s < M* w ^ h p 
(or equivalently such that h Ps < M s w and Ap g+1 > Mffi 1 ). 

(5) We calculate f Pr , r= 1, . . ., ra, using formula (2.12). 
Suppose now that weW with & w disjoint. We define Effl by 

E&& = &' (2.19) 

where &' is the solution of the minimization problem (2.7) for the given ,W . This definition induces 
the definition of Eaaj for simple networks with quadratic cost functions. 

Note that (2.7) simply states that ,W =Ep^ satisfies (1*13) for the pair w. Thus, if EffiP 
= <!F, J^ satisfies (1.13) for the pair w, i.e., Effl satisfies condition V of Theorem (2.2). Condition 
2' of the same theorem follows immediately from the continuity of the functions involved in (2.12). 
Finally, conditions 3'. 4' are also satisfied since SP ' = E {[V, F minimizes C(J^' ) over the set 3? ,/ 

By Theorems (2.1), (2.2) it follows that E (isj induces an algorithm for a simple network with 
quadratic cost functions. 

We should emphasize here that the effectiveness of the solution of the minimization problem 
(2.7) is due essentially to the assumption of the disjointness of &* w and of the quadraticity of 
{c a (fa):ae / }. 

We now proceed to extend the definition of Effl in the case for which (P w is not disjoint. Note 
that the minimization problem (2.7) is meaningful even in this case. The solution of this problem 
would provide a natural extension of the definition of Effl to cases where £P W is not disjoint. Un- 
fortunately an effective solution of this problem does not seem possible. Thus we devise the follow- 
ing kind of extension: 

We choose to calculate 

:F' = Ep^ 

by following the steps (1), (2), (3), (4). and (5) described above. 

This is clearly possible since this procedure, though motivated for disjoint lP w , does not de- 
pend in its definition on the assumption of disjointness. In this way we retain the simplicity of the 
calculations. Let us show that the condition 1 ' of Theorem (2.2) is satisfied. Suppose that 

M w — h p 

&'=E W & = & for some ;¥*%. From (2A2),f Pr =fi = —z -, r= 1, . . ., s, and/ Pr =/;=0, 

*gp r 

r=s+l,. . ., m. Then, recalling (2.9) and (2. 10), fip r = M W9 r=l 9 . . .,s,/%> r >M w ,r=s+l, . . .,/n. 

Now observe that yu />r is the (real) marginal cost along the path p r . Hence the equilibrium equations 
(1.13) are satisfied for w and the proof of condition 1' is complete. Furthermore, we can prove, as 
before, that the condition 2' of Theorem (2.2) remains valid. On the other hand, the motivation 
which was present in the case of disjoint ^ w and which was justified by the proof of the validity of 
conditions 3', 4'. of Theorem (2.2) is not present any more. Thus, it is not obvious that conditions 
3', 4' are still satisfied, and we have to go through the following lengthy calculation in order to 
check whether they are valid. 

Let us set AJ*~= ^' — :W . The change of the total cost is given by 

C(F')-C(¥) = J t [ga(f„+£~f a ) 2 + h a (f a + l~f a )- ga fz-hja}= 2 gAfl 

+ 2 [2gJ ri + /,„]!/„= 2 A/p X [2gja+h„]8 ap + X ft»4/2 

(ie/ pe&w aeJif ae J 

= 2 A/„[2 ft ,/„ + /,„| l ]>>Af« 
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where use has been made of (2.9). Thus 



By (2.11) 



C(*")-C(^)=-2 2 «pA/J+ X Vh[2&/;+Ap] + 2 -" A/ ''- 

pe^ 5 pe^ n« / 



8 /« IK 

2 A/„,.[2 ft ,/;, + M = 1/,, 2 A/,,- 2 (~ WW § ¥, 2 A/ Pr =0. 

pt &> r=l r=s+l r=l 



(2.20) 



Henc 






Recalling the Definition (2.8) of g p < and (1.4), 



C(P')-C(&) ^ 2 g a 2 BaAfp " 2 S 8 ^A/| 



pe^«- 



(2.21) 



We want to study now the sign of C{^F') — Cffi). It is convenient to consider links of a special 
type. A link a will be called simple, doable, or total with respect to w depending upon whether a is 
contained in precisely one. two or all of the paths of & w . Note that if !? w is disjoint, then all links 
contained in £P W are simple. 

Suppose now that a is a simple link contained in a path pe^ w . Then, 



g„\\ £ 8„ p A/ p T-2 2 8„A/|] = -g„A/Jg0 



(2.22) 



with equality holding if and only if kf p = 0. 

Let a be a double link contained in two paths p, qe& w . Then, 



go 



2 8 ap Af p -2 £ 8„ p A/|[=ga{[A/ p + A/ g ] 2 -2[A/|-fA/2]}=-^(A^-AA)^0 (2.23) 

pe^V -1 pe^V 



with equality holding if and only if Af p = Af (j . 

Suppose finally, that a is a total link with respect to w. Then, 



g n l\ 2 8 OJ) A/ p T-2 2 K^A=~2ga 2 A/|£0 

^L peA- -" pe^«- ' pe^V 



(2.24) 



with equality holding if and only if Af p = 0, pe^ w . 

Thus (2.21) implies that, if all links contained in £P W are simple and/or double and/or total 
links, then 



C{^')-C(^ )^0. 



(2.25) 



Two paths of $P W will be called directly connected if they share at least one double link (with 
respect to w); will be called connected if they can be joined by a finite sequence of direct connec- 
tions. We now introduce the following definition. 

DEFINITION (2.2): The set & w will be called almost disjoint if the following conditions are 
satisfied. 

(1) All links contained in £P W are simple, double, or total with respect to w. 

(2) Any two paths of J^ r , which consist exclusively of double links and are connected only to 
paths consisting exclusively of double links, must be connected to each other. 
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If 3P W is almost disjoint for all weW\ the 1 transportation network will be called a/most simple. 
In particular a simple network is also almost simple. 

Assume now that & w is almost disjoint. Then (2.25) is satisfied. In addition, we claim thai 
equality can hold only if A<^"=0. In fact, suppose that equality holds in (2.25). From (2.22). (2.24) 
it follows that 

A/ p =0 

in the case where pe^ w contains at least one simple and/or total link. Let p, qe^ w be two paths 
which are directly connected. On account of (2.23). 



A/ P = A/* 



(2.26) 



It is clear that then the validity of (2.26) extends to the case where p and q are connected. Thus 
4/p^O even for all pe£P w which are connected to a path containing a simple and/or total link. In 
,>rder to complete the proof it remains to consider the set ^/, of paths in {P w consisting exclusively 
of double links and not connected to any path containing simple and/or total links. On account 
of the definition of almost disjointness, if p, qe^' w , then p and q are connected in which case 



Recall the conservation equation 



a/„=a/;, = a/. 






(2.27) 



(2.28) 



Using (2.27), (2.28) we deduce that A/= 0, i.e., (2.25) holds as equality only if A.^^0. 

Combining the above results with the Theorems (2.1), (2.2) we reach the following conclusion. 
THEOREM (2.4): The equilibration operator Edsj induces an algorithm which solves problem 
P\ for almost simple networks, with quadratic cost functions. 

Actually, this result is the best possible in the sense that there exist networks with triple 
links for which the operator Edsj does not induce an algorithm which solves problem P\ for arbi- 
trary initial ;W m . An example of such a network is given in appendix 2. Nevertheless, recall that 
EjV satisfies conditions 1', 2' of Theorem (2.2) for arbitrary 3P W . From Theorem (2.3) it follows that 
if {E' ( \ S j:F m ) converges then it will converge to the solution J^, of problem P u In consequence, 
it is worthwhile to try an application of E dsj even for networks which are not almost simple. 

In order for an algorithm to be appropriate for application in practice, it is not enough that 
it converges; it must converge rapidly. We have the following evidence about rapid convergence 
of the algorithm induced by Edsj in the case of almost simple networks. We consider the test net- 
work of figure 1 (with 60 paths) which is almost simple, but is not simple. 
We have developed a computer program which solves problem Pi for this 
network by application of Edsj- We have calculated the solution for a very 
wide range of demands and choice of the initial distribution. We have 
observed extremely rapid convergence. More specifically if g, h, were 
the average values of ga, ha, aeJ£, and d is the average demand in the 

c - — 

network, we chose -entering in Definition (1.2) equal to 2gd+h, and after 

5 iterations the flow pattern satisfied eqs. (1.13) modulo 10~ 5 . Furthermore, 
we have treated the same network algebraically and have shown that 




Figure 1. 



||£^ (0) -£V^ (0) || 



decreases with the speed of a geometric progression with a ratio less than 1. 

With these observations we conclude the discussion of the equilibration operator Edsj for 
transportation networks with quadratic cost functions. 
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b. The Equilibration Operator E n dsj for the Quadratic Model 

In this section we introduce an equilibration operator E n dsj which does induce an algorithm 
for the solution of problem Pi for arbitrary transportation networks with quadratic cost functions. 

Fix some weW. The motivation of the introduction of E T $ sj is similar to that of E*W. Namely, 
Endsj w jV} b e selected to be a minimization operator of C(J^) but over a set less broad than the 
set J3?w defined earlier. 

Let <jFeiF. We define &' = E r $ sj ^F by the minimization problem 

minimize C(J^'), over ^'e&w where Jf^, to be defined, is a subset of 2£$. containing J*\ 

(2.29) 

Note that the minimization problem (2.29) is very similar to (2.7). Such a definition of Effii 
guarantees automatically that condition 3' of Theorem (2.2) is satisfied. The main difficulty is 
that we must select the set iz ifin such a way so that conditions 1', 2', 4' of Theorem (2.2) are 
also satisfied while, at the same time, the solution of (2.29) can be obtained in an elegant way. 7 

The marginal cost on a path teSP w corresponding to the flow pattern !F is given by \x Vr {3F) 
of (2.9) with p r =t. 

Let p, q €& w be defined by 

(2.30) 
/V^~) = max {iLt(&)it& W9 ft > 0} 9 . 



We define the set iff by 



£Z = {,t'e£:f t =f t , unless t = p % q}. (2.31) 



We now define E™W J* = J^'as the solution of the minimization problem (2.29) over & w as 
selected above. 

The calculation of &' =Efi, d8J ^ amounts to the calculation of the two new flows f p ,f q . 
Let us define 

4=2 %?<» «f s 2 *»*« < 2 - 321 

aeJf a*Jf 

where the incidence symbols 8% p , 8f, q have been defined by (1.19). 

The minimization problem (2.29) leads to the following equilibrium conditions analogous to 
(2.11): 

lip(P)+2&(f p -f p ) g fju q (^) +2g$ (f' q -f q \ 

f' P +n=h+U (2.33) 

with the understanding that, if (2.33)i holds as a strict inequality, then the second one of (2.33) ;} 
must hold as an equality, and conversely, if both (2.33) 3 hold as strict inequalities, then (2.33)i 
must hold as an equality. The solution of (2.33) is given by 

Jp jp^ 



2(gJ+g§) 



7 If we selected S w — £*,■> we would be led to the minimization problem (2.5) v In this case we know that (2.29) does not have an elegant solution unless ^V is 
disjoint. This observation emphasizes the fact that the selection of an appropriate 3?&\% not easy. 

8 If n attains its minimum for more than one path, then select p as any of those paths. 

9 If ix attains its maximum for more than one path, then select q as any of those paths. 
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Jq J " 2($+g$) ' (2.34) 

17 "- 2(«J+«f) ' 

or by 

f P =f P +f g , 

Jl = 0, (2.35) 

11/9 2(g«+^) 

Thus, the selection ofiFifhas led to a minimization problem (2.29) whose solution can be 
calculated efficiently through (2.34) or (2.35). 

We now proceed to prove that E^ satisfies the conditions l'-4' of Theorem (2.2). From (2.33) 
it is clear that E" { ^:¥ = ^ if and only if fi p (^) = fJL q (^) (recall that /,>()). 

The paths /;, q have been selected so that for any te0> w with/* > 0, 

Thus, if E^ 8 ^ = $F , then all paths £e,^, r with// > have the same marginal cost and the equilibrium 
conditions are satisfied for w. Hence Ej' r (h J satisfies condition 1' of Theorem (2.2). Condition 2' is 
obviously satisfied. 

Let us now calculate AC = C(jF) — C{E\\! h >^F). Obviously AC is given by 

AC = ^ 8% \ goft + h qf a - gcfa 2 - hja ] + ]£ 8g g [gj% + ft /a ~ g«/a 2 ~ ^ /a 1 

= ^ 8*p[W« + /,„ - g„(f a -/;, ) ] (/;, -/; ) + x s/; " t 2 ^- + /; " -«*& ~f« > ] (£ -/«)• 

Note that in the first sum/ a — fa—fp~ fp-> while in the second sum / a — /« = fq~fq- Hence, recall- 
ing also the definition (2.32) of g#, gf, we get 

AC= (fp-fp')li$-gp 1 (fp-fp') 2 + (f q -fq)l4~^(f q -f q ) 2 (2.36) 

where 

M* s 2 ^Rft/a + VU (2.37) 

K S E 8g,[2^/„ + A a ]. (2.38) 

By (2.33),,/, -J? — (/;,-/;)• Hence 

ac = [ fit+^w; -/, ) - k - «c < /,; -/, » ] ifp -/; ) 

= [/i|+2g|(f;-/,)-M€-2«s(/« -/«)](/* -/;)+(^+«f)(/p-/;) 2 - (2.39) 

But /a# — ju.f ; ' = ju. /; — /x (/ . Then, using (2.33), and the fact that/ p — f' p S 0, we obtain 

AC ^(^ + g§) (/,-/,;)* & (2.40) 
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with equality holding only if fp=f P , i.e., if E^^=^. Hence conditions 3', 4' of Theorem (2.2) 
are satisfied. Combining the above results we reach the following conclusion: 

THEOREM (2.5): The equilibration operator E nds J induces an algorithm which solves problem P x 
for an arbitrary transportation network with quadratic cost functions. 

We wish to compare the algorithms induced by the equilibration operators Edsj, E n dsj and to : 
point out their corresponding advantages. We prefer to postpone this comparison until Edsj, E m i S j 
have been extended to cover the case of nonquadratic cost functions. This extension is the subject 
of the next section. 

2.3. The General Model 

In this section, we will extend the definitions of the equilibration operators Edsj* Endsj to the 
case of a general transportation network 2T = (^, £d , c €). We will assume that c a is twice continuously 
differentiable for all ae^f. 

Note that the minimization problems (2.7), (2.29), through which the operators Eff, E r ^ sj have 
been introduced for the quadratic model, are well set also for the general model. Thus it appears 
that the proper extension of Effl. E%* s i would be obtained through the same minimization prob- 
lems set for general cost functions. However, a review of the theory described in section 2.2 indi- 
cates that the simplification induced by the assumption of quadratic cost functions lies in the fact 
that for such cost functions the equilibrium conditions (2.11), (2.33) corresponding to (2.7), (2.29), 
respectively, are linear. This fact permits a very effective and elegant solution for both minimization 
problems (2.7), (2.29). 

In the case of general convex cost functions, the equilibrium conditions are in general non- 
linear and hence a simple solution of problems (2.7), (2.29) is no longer possible. In order to devise 
a working extension of the definition of Edsj, Endsj for the general model we use the following 
considerations: 

By Taylor's theorem, if f' a \s close to f a , 

C„(f' n ) = Ufa) +C' a (f a ) (/'„-/„) +ic'a(fa) (/W«)*- ' (2.41) 

In particular, for quadratic c„(/„) = g l{ f„ 2 + h,f„. 

Ca(fa) = C„(f„)+ (2gaf„ + h a )\f'„-f„) + gaifi-fa)*. (2.42) 

Comparing (2.41) with (2.42) we observe that the quadratic cost function ga(fa)f'a + h a (fa)fa which 
approximates the general cost function c u (fa) in the neighborhood of/,, has coefficients 

g n (f„)=Wa(fa), (2.43) 

ha(fa)=c'a{f ll )-cl(fa)fa. (2.44) 

Using as a motivation this observation, given SFe^E we construct the collection of quadratic 
cost functions 

«- Hc'a (fa)fa i + Waif a) ~ cZ(f„) falfcacS?} . (2.45) 

If 

We now define E d t p (resp. EffW) by identifying Ew>3F (resp. E n $ 8 ^) for the transportation net- 
work 5~= (^, 0, #) with E$& (resp. E'tf*^) for the transportation network ^ = (#, 2), %?), 
the latter being well defined since *§V consists of quadratic cost functions. In other words, 
&' = Effi3F (resp. &' = Eff"^) will be calculated through (2.12), (2.13) (resp. (2.34) or (2.35)) 
where 

gp(^)^\^8apd'jj a ), (2.46) 
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M-^)= ^^apC'aifa), (2.47) 

hp(&)=n P (& r )-2g p (&')f p < (2.48) 

tiWml^fycTJja). (2.49) 

Hence the definition of Edsj< E n dsj has been extended to transportation networks with general 
convex cost functions. 

The numerical application of these operators to a given ,Fe& follows precisely the same 
pattern as in the case of the quadratic model, the only difference being that the coefficients g p , 
h { ,, fi p . ,gp* are not constants any more but should be calculated at every step. 

It remains to examine under what conditions Ed S j, E n dsj do induce algorithms for the solution 
of problem l\. As in the case of the quadratic model, conditions 1', 2' of Theorem (2.2) are satisfied. 
It is intuitively expected that 3', 4' of the same theorem are more apt to be valid if the approxima- 
tion of % by Sjf is faithful, that is if c a (fa) is sufficiently close to a quadratic function, aeJ/£ . This 
idea has been verified for simple networks in [10], where we give sufficient conditions so that 
Effl satisfies 3', 4'. We will go here through a similar analysis for the operator E','/ 1 ^ in the case of 
a general network. 

Fix ire//' and let ./c/. We set &' = Etf 8 ^. The change of the total cost 

M = C(W)-C(W') (2. so 

is given by 

AC=2 8|p[c„(£)-C«(/2)]+ X $& q [Ca(fa)-Ca(f«)] (2.51) 

<ie/ <i€ / 

where p, qe./\ r are the paths which are "equilibrated." Applying Taylor's theorem, 

M: = ~ S ^U ( 'n(fa) (fa -/*) + « Ci(f„)(fl -fa) 2 ] 

at y z 

- 2 ««</ W (/« ) Ha -fa ) + g Ca(fa ) (fa ~fa ) 2 ] 

<ie y 

with f„ between f a and /',. Note that in the first sum fa—fa = fp~fp while in the second sum 

fa ~~fa ~~Jq~Jq- 



Let us set 



K' = rfW^^ 8?;, c,; (/„ ) , (2.52) 

ae / 

n = M| ( .^ ) ■ £ 8 V« (/« > ' (2.53) 

g?,=m&)=l 2 S&fi'aifa). (2.54) 

tt = g${&)=\ 2 KCaifa), (2.55) 



and recall that / ; , —j], = — (f q —f' q ) . Then 
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= \j4 + 2g$w P -f„) - ixij - 2§f (^ -./;,) ] (/ P -/; » 

+ [2g> - £J + 2g* - gP] (f v -/; ) 2 (2.56) 

where #$, #§ stand for gJ(«!F), g^i.W) as denned by (2.49). On account of the definition of E$W 9 
we liave the equilibration condition (2.33)i written in the form 

pj + 2g« (f P -f;)-fiP q -2gP (f q -f' q ) ^0, (2.57) 

where use has been made of the obvious equation 

Combining (2.56) with (2.57) and the fact thatfp— f p iOwe obtain 

&C^ [2g$-g$+2g%-g q >](f;-f p ) 2 . (2.58) 

From this last result we deduce that conditions 3', 4' of Theorem (2.2) are met if 

2g9(^)-g9(#)H-2^(^)-^(#) >0 (2.59) 

for all possible «^*, j^. 

A sufficient condition for the validity of (2.59) is that 

c" n {f„) >\c"„(fa) (2.60) 

for all possible f a ,fa, aeJ^. Suppose that we know that/ rt ,/^, are limited in some interval /„. Then 
f a will also be limited within I„. Condition (2.60) is obviously satisfied if 

minc';,(/)>-!maxc;(/). (2.61) 

'a & l a 

Thus, if (2.61) is valid for all ae^, then E#W satisfies conditions 1', 2', 3', 4' of Theorem (2.2) for 
any weW. Using Theorems (2.2) and (2.1) we reach the following conclusion: 

THEOREM (2.6): Let 2T= (& y 3) y %) be a (general) transportation network. If condition (2.61) 
is satisfied for all aeJ^, then the equilibration operator E n ^ S j induces an algorithm for the solution 
of problem P 2 [&~]. 

Actually condition (2.61) states that the oscillation of the function c a f (f a ) on I a is not very large, 
or in other words, that c (l (f n ) is sufficiently close to some quadratic function on /<*. Thus Theorem 
(2.6) is in accordance with the intuitive idea cited before. 

In order to put the assertion of Theorem (2.6) into practical use we have to find intervals I a 
with the propertv stated above. 

Using the feasibility condition (1.1) we conclude that we may take 

/«= [</>„, *«] (2.62) 

where 

wS r a weW n 

114 



Here W a stands for the set of all weW which are connected by at least one path p containing the 
link a, and^, stands for the set of a\\weW (l with respect to which a is total (if W'a = J, then (f> a = 0). 

Note that (2.61) is more apt to be satisfied if I a is "small." Actually the /„, as given by (2.63), 
are the best possible (i.e., the smallest possible) if we are to expect a convergent sequence 
{E^asj^^} f° r every ^ (()) e^. In practice, though, we are merely interested in knowing whether 
{Endsj^^} converges for a specific J^ (0) , namely the one selected as the starting point. If we re- 
strict ourselves to this problem, then it is possible, at least for special types of networks, to obtain 
I a which are proper subsets of the I a given by (2.63) and hence are preferable. 

Concluding this section, we want to emphasize that a very large subset of the set of cost func- 
tions which satisfy conditions 1, 2, 3, 4, of section 1.3 do satisfy also (2.61), and hence the operator 
Endsj induces an algorithm for the solution of problem Pi for a very wide class of transportation 
networks. We also want to emphasize that it is worthwhile to try to solve problem P\ by means of 
one of the operators Edsj or Endsj even if conditions 3', 4' of Theorem (2.2) are not met. In fact we 
have shown that the above operators always satisfy conditions 1, 2 of Theorem (2.1) and hence, 
according to Theorem (2.3), if {E n ^ {{)) } converges, then it will converge to a solution of problem P\ {) . 

2.4. Comparison Between Edsj and E n dsj 

From the theoretical point of view E m t S j is superior to Edsj since it can be used for the solution 
of problem P\ even in the case of not (almost) simple networks. From the practical point of view, 
though, Edsj also has some advantages. In fact Ep equilibrates the whole set of paths 8P W but does 
so perfectly only if P? w is disjoint. On the other hand E^ equilibrates only the two "most un- 
balanced" paths but does so perfectly in all cases. It is clear that in the case of a network which 
is not almost simple we must apply E m tsj. 

In the case of an almost simple network we advise the application of Endsj if most of the links 
are simple and @> w contains many (more than two) paths for at least one weW. On the other hand, 
we advise the application of Ed.y in the case of a network in which 0* w contains few (of the order 
of two) paths for all weW and there are relatively numerous double and/or total links in the network. 

3. Appendix I. A Network With a Nontrivial R[W] 

We consider the network of figure 2 with links 

a\= Ui, * 2 ), a 2 = (*2, * 3 ), a 3 = (x u x 4 ), a 4 = (xt, Jt 3 ), 

a 5 = (*3, *5), a 6 = (x 5 , x 7 ), a 7 = (* 3 , x 6 ), 08 = (*e, x 7 ) , 
and paths 

pi= (a u a 2 , a 5 , as), IH= (a u a 2 , a 7 , a 8 ) , 

p 3 = (a 3 , a 4 . a 7 , a 8 ), P4 = (a 3 , "4. a 5 , a%). 

Let d be the demand associated with the pair (xi, xt). For this network 
we have 

J Pi +JP2 ~^~fp3 ~^~fp 4 ~ d, 

J Pi ~^JPz ~Jai , fp x +fp 4 ~fa 5 , 

J Pi ' J Pi = Ja 2 , J Pi ~T~fp 4 =z Ja 6 , 

fp3~^fp4~fa a > fvt+fpz—fai, 

Figure 2. 




u Such examples have been constructed. 

115 



/ Pf ^0,5 = 1, . . ., 4. (Li) 

The feasibility conditions on 3F read: 

fa g ^ 0,5=1, . . .,8, 

7 «i = /«2' /«» = fa 4 '> fas == /a 6 i /a 7 = Ja H -> \*-^) 

Ja x ~rfa 3 = /« 5 +/o 7 = «• 

Suppose now that J^ is given such that (1.2) are satisfied. 

Assume first that at least one of the f a is 0. Then it can be shown that there exists a unique 
solution of (1.1), i.e., in this case, R[W] contains a unique element. 

On the other hand, if ZF is such that f a > 0, 5 = 1, . . .,8, then it can be shown that the 
solutions of (1.1) form a one parameter family, i.e. /?[iF] is a convex subset of a one dimensional 
vector space. 

4. Appendix II. An Example of a Network Which is not Almost Simple and for Which 

Edsj Fails 

In this appendix we present an example which shows that the assertion of Theorem (2.4) is 
the best possible. To be precise, we consider the network of figure 3 with links 



ai= (xi,x 5 ) 9 a 2 = (x u x 2 ),a 3 = (^2,^5), a 4 = (^2,^3), 

«5= (#3, #5), a 6 = (x 2 ,X4),a 7 = (X4,X 5 ), 

and paths 

pi = a u P2= («2, «4, as), P3= («2, a 3 ), 

P4= (a 2 ,a 6 ,a 7 ). 




Figure 3. 



Let d be the demand associated with the pair (x\, x 5 ) large enough so that all paths operate 
at a nonzero level. We assume that the cost functions are quadratic of the form 

Ca(fa) =gafi + h a fa, 

where ga x = ga 3 = 2g, g tt2 = 8#, g tt4 = ga 5 — #a 6 — ga 7 — g, 

h (h = 3h, h (l2 = ha 4 = ha 5 = h (H =ha 7 =h,ha 3 — 2h, 

with arbitrary g > 0, h ^ 0. 

Note that the link a 2 is triple and hence this network is not almost simple. Thus Theorem 
(2.4) does not guarantee that {E'^j J^ (0) } converges. In fact we will prove the following interesting 
result. The sequence {E% sj J^°>} converges only if the initially chosen J^ (0) satisfies /^ ) = /i Pl where 
<jFi is a solution of problem P\. 
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To prove this result select J^eiF. We want to calculate Af Ih =fj, 1 — f f > r It will turn out below 
that h Pl is the smallest of the /i, v From (2.10), (2.12), 

A/,„ = -^. (IL1) 

Proceeding to the calculation and using (2.8), (2.9), (2.10) and the obvious relations 

J<t\ = /Pu fa* —Jin *JPz ~^JP4-> Ja 3 = JP^ Ja 4 = fa 5 = Jp 2 ,Ja* = fa 7 = fp 4 , 

we end up with 

gpi = 2g, gp 2 = gPz = gp 4 = 1%, 

h Pl = 3h, 
h P2 =l6g(f P4 +f P3 )+3h, 
h P3 =l6g(f P2 +f P4 )+3h, 

h P4 =l6g(f P3 +f P2 )+3h. 

Whence, from (2.13), 

Ugd+6h-8gf Pl 
2 
where use has been made of 

Substituting into (II. 1) we obtain 

A/p= iM^. (IL2) 

Suppose now that we apply E (lsj on &\ Let h,^' = E d ^' —&'. From (II.2), 

13 ( i- 16,/;, _ 13 ( i-16(/„ + A/„) _ 13^-16/^ _ __ 

4//>, = g ~~ g ~~ 8 ^A/Pi — A/Pi mjpi— a JPi- 

Thus, {E>l Uj tF i()) } will cycle unless <F (()) has been chosen in such a way that A/JJ^O. 11 

5. Appendix III. Calculation of the Smallest Number Modulo Which a Flow Pattern 

Satisfies the Equations of Equilibrium 

In this appendix we present a procedure which enables us to calculate the smallest number 
e module which a given flow pattern J^ satisfies the equilibrium equations (1.13). 12 

Fix weW, and let/i, . . ., f m be the flows on the paths of @* w and c{, . . ., c' m the correspond- 
ing marginal costs. We want to find the minimum of the values of e^ such that for any p, q in 
{1, . . ., m} , if fp 5* de u >, then either 

I c p~ c q I < ce wld holds, or both of c' p < c' q + ce w /d an\df Q < de w hold. 



1 Note that if A,fffl=Q we can apply Theorem (2.4), since a 2 becomes then total for the network emerging by the omission of the path p\. 

2 This method has been kindly communicated to us by Alan Goldman. 
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This last disjunction is equivalent to an exclusive disjunction: either | c' p — c' q \ < ce w /d holds, 
or f q < de w and c' p — c' q < ce w /d ^ | c' p — c' q | hold, with the last condition equivalent to the conjunction 
of c' p < c' q and ce w /d ^ c' q — c p . 

Number so that 0=/o^/i ^ . . . ^f m . We ask whether it is possible to choose e w in the 
interval (ft/d, ft+i/d]. For such an e w , the requirement is that if p > t, then for each q either 
e w > (d/c) \c' p — c'q | , or q ^ t and c' p < c' q and e w ^ (d/c) (c q — c p ). That is, if p > t, then (a) q > t 
implies e w > (d/c)\c' p — c q \, while (b) q^t and c' p ^ c' q imply the same conclusion. Since this is 
to hold for all p > t, we see that (a) is equivalent to 

€ w > (d/c) max{ \c' p — c' q \ :p>t, q>t}= (d/c) [max p>t c' p — min 9>/ c^] , (HIT) 

while (b) is equivalent to 

€^> (d/c) mdLx{c' p -c q :q^t< p, c' p ^ c q }= (d/c) max{max p >^- mm Q ^ t c' q , 0}. (III.2) 

Combining (III. 1) and (III. 2), and noting that the lower bound for e w in (III.l) is nonnegative, we get 

e ir > (d/c) [max„>^,— mm q> tc q , mdLX p>( c' p — m\nq <t c'q] = (d/c) [max p>( c' p — c^ in ] (III. 3) 

where c^in^mm {c\, . . ., c' m ) . But (III. 3) can hold, for an e w in interval (f t ld,f t + i/d], if and 
only if 

ft + i/d>(d/c) [max p>t c' p -c^ n ]. (III. 4) 

Thus the infimum of the allowable e ^-values is f t /d for the smallest t such that (III. 4) holds. 
We set 

e = min e w 

we.W 

and obviously this e is the smallest number modulo which (1.13) are satisfied by 3F . 
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